Static and dynamic properties of crystalline phases of two-dimensional electrons 
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We study the cohesive energy and elastic properties as well as normal modes of the Wigner and 
bubble crystals of the two-dimensional electron system (2DES) in higher Landau levels. Using a 
simple Hartree-Fock approach, we show that the shear moduli (cqq's) of these electronic crystals 
show a non-monotonic behavior as a function of the partial filling factor u* at any given Landau 
level, with cee increasing for small values of u* , before reaching a maximum at some intermediate 
filling factor u^, and monotonically decreasing for v* > v* n . We also go beyond previous treatments, 
and study how the phase diagram and elastic properties of electron solids are changed by the effects 
of screening by electrons in lower Landau levels, and by a finite thickness of the experimental sample. 
The implications of these results on microwave resonance experiments are briefly discussed. 
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I. INTRODUCTION 



There has been much interest recently, both theoreti- 
cally and experimentally, in the quantum phases of the 
two-dimensional electron system in higher Landau levels 
(LLs). Hartree-Fock (HF) studies^ have shown that, for 
small partial filling factors v* = v—2n 0.1 — 0.2 (with v 
the total filling factor and n the Landau level index), the 
electrons form a triangular Wigner crystal (WC), while 
for v* close to 1/2 the ground state of the system is 
a unidirectional charge density wave or "stripe" stated 
Between these two regions, HFi&a and density matrix 
renormalization group (DMRG) methods, 5 as well as ex- 
act diagonalization on small systems^ have suggested the 
existence of a new kind of crystal structure, the "bub- 
ble crystal" (BC), with more than one electron per lat- 
tice cell. While HF predicts transitions between bubble 
phases with increasing number M of electrons per bub- 
ble, with M going up to M = n + 1 electrons in LL n, in 
DMRG the last bubble phase is absent, and a transition 
takes place directly from the BC with M = n electrons 
per bubble to the stripe state. It is noteworthy that the 
difference in cohesive energy between all these different 
phases remains very small, of order 0.01 — 0.03e 2 /e£ (with 
e the electronic charge, e the dielectric constant of the 
host semiconductor and £ the magnetic length). 

On the experimental side, evidence for the existence of 
these different phases is found mainly through transport 
experiments. Early DC measurements^ have shown 
strongly anisotropic transport around half-filling, which 
was interpreted as evidence for a stripe state. Crystal 
phases of electrons, which are pinned by quenched dis- 
order, should all be insulating in the DC regime. One 
has therefore to resort to other experimental techniques 
to resolve the transitions between these different crys- 
tal phases (e.g. WC to BC). One particular such tech- 
nique consists in measuring the microwave response of the 
2DES, which would give different resonant behavior for 
different electronic crystalline phases. While early mi- 
crowave conductivity experiments^*^ gave considerable 
support to the existence of a WC around integer filling 



v = 1,2,3 and 4 for v* < 1/2, it was not until recently 
that coexistence between two phases with two distinct 
resonant peaks was observed by Lewis et ali^. It is our 
aim in this paper to discuss in some detail the experi- 
mental results of this last reference, in light of recent ad- 
vances in resonance pinning theories of two-dimensional 
electronic solids JSil/L In order to do so, we will need to 
develop a microscopic picture of the 2DES as a system 
of interacting guiding centers. Indeed, after projection of 
the electronic density on the uppermost Landau level the 
problem at hand reduces to one of interacting guiding 
centers (of quenched kinetic energy) with an appropri- 
ately defined effective interaction potential that includes 
quantum effects at the Hartree-Fock level. The knowl- 
edge of the effective interaction potential between guiding 
centers will allow us to find the elastic moduli and the 
normal eigenmodes of the Wigner and bubble crystals. 
In agreement with the results of Maki and Zotosi^ for 
the WC in the lowest LL, we find that the shear moduli 
((Ws) of the Wigner and bubble crystals show a non- 
monotonic behavior as a function of the partial filling 
factor v* at any given Landau level, with cqq increasing 
for small values of v* , before reaching a maximum at 
some intermediate filling factor v* nl and then monoton- 
ically decreasing for v* > v* n . These results will allow 
us to attempt a qualitative analysis of recent microwave 
experiments by Lewis et alm^. We find that, while the v* 
behavior of the first resonance peak observed in this last 
reference is in good qualitative agreement with existing 
theoretical predictions for the pinned Wigner crystal^ 
the second resonance peak has a behavior as a function 
of v* which is quite different from what a linearized so- 
lution of self-consistent equations for the microwave re- 
sponse theories predict for the two-electron bubble state. 
Below, we shall argue that an adequate description of the 
resonance peak of bubble crystals may require the full nu- 
merical solution of the self-consistent response equations, 
as done in Ref. 24. 

This paper is organized as follows. In Sec. |ll|we review 
the physics of the Wigner crystal and bubble states, and 
introduce some notation. In Sec. IIIII we find the shear 
moduli of the Wigner and bubble crystals in LLs n = 2 
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and n = 3. In Sec. IIVI we calculate the normal modes of 
the Wigner crystal and of the two-electron bubble solid 
in the Landau level n = 2. Then, in Sec. Elwe investigate 
the effect of screening by lower LLs and of finite sample 
thicknesses on the cohesive energies and shear moduli of 
the Wigner and bubble crystals in LLs n — 2 and n = 3. 
In Sec EH we discuss the implications of our results for 
the shear moduli on the microwave response of electronic 
crystals in light of the recent microwave conductivity ex- 
periments of Lewis et alM> Finally, Sec. IVIII contains a 
summary of our results along with our conclusions. 



II. WIGNER CRYSTAL AND BUBBLE PHASES 
IN HIGHER LANDAU LEVELS: 
HARTREE-FOCK APPROACH 

We shall start from the expression of the Hartree-Fock 
Hamiltonian of the n-th partially filled Landau level, 
which is given by (throughout this paper, we use J as a 

shorthand for / 75^) 



H H F = \fv HF (c i )\p(c i )\\ 



(1) 



where Vhf is the Hartree-Fock interaction, and p(q) is 
the projection of the electronic density n(q) onto the up- 
permost LL, and is given by 

P(q) = e -<^ 2 /4 Ln (^2 /2) - W 

For a bubble crystal with M electrons per bubble, we 
shall approximate the electronic density by 



M-i 



n 



T) = E E K m (r-R,)| 5 



(3) 



i m— 



where ip n ,m(r) is the noninteracting wavefunction of an- 
gular momentum m and Landau level index n, and where 
the summation extends over all the bubbles located at the 
lattice sites of a triangular Bravais lattice and over all 
the electrons within each bubble. The Hartree-Fock in- 
teraction potential in Eq. Q consists of the sum of a 
Hartree and Fock parts, which are given, in the Landau 
gauge A = Bxy (A being the vector potential), by&i 

Vff(q) = — e-^ 2 / 2 [L n (^ 2 /2)] 2 , (4a) 



eq 



V>(q) = -(2^ 2 ) f y ff (q')e- lqx ^ 



(4b) 



where L n (x) is the n-th Laguerre polynomial. 

The calculation of the cohesive energy of the Wigner 
crystal and bubble phases proceeds in a standard way as 
follows. From Eqs. 10 and (0, it is easy to see that we 
can write the projected density p(q) in the form 



p(q) = p(q) ' 



(5) 



where the sum runs through all vectors of the trian- 
gular Bravais lattice, and where p(q) is given by 



p(q) 



n(q) 



-9 2p /4 iri ( (? 2^2 / /2)' 



(6) 



with n(q) the Fourier transform of the density at a given 
bubble, 



M-l 

(q) = XI / dv Wmn{v) 



V iq ' r . 



(7) 



Inserting the decomposition © into the expression of 
the cohesive energy, Eq. and making use of Poisson's 
summation formula (here A c — 2iri 2 M/v* is the area 
of the unit lattice cell and the Q's are reciprocal lattice 
vectors) 



£e-^ = ^5> (q -Q), 



(8) 



Q 



we finally obtain the following expression for the cohesive 
energy per particle E coh = (H HF )/N : 

Eco h = [^(Q)(i-^q.o) + ^(Q)]Ip(Q)I 2 , (9) 



Q 
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FIG. 1: (Color online) Cohesive energy of the Wigner crystal 
(solid line), M = 2 (dotted line), M = 3 (thin dash-dotted 
line) and M = 4 bubble solids (thick dash-dotted line in the 
lower panel). The cohesive energy of the stripe state is indi- 
cated by a dashed line, and all cohesive energies are in units 
of e 2 /el. The upper panel is for n = 2, and the lower panel is 
for n = 3. 
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where charge neutrality requires that we exclude the di- 
verging Q = Hartree term from the summation. In 
Fig. ^ we show the cohesive energies of the different 
phases (Wigner crystal, bubble solids and stripes) in the 
Landau levels n — 2 and n — 3. The values of the cohe- 
sive energies shown, and of the partial filling factor v* , 
at which transitions between the different phases occur 
in both Landau levels, are in excellent agreement with 
recent Hartree-Fock calculations by Cote et al~ and by 
Goerbig et alX As it can be seen, the energy difference 
between the WC, BC and stripe states is extremely small, 
of order 0.02e 2 /e£, which is within the range of the mag- 
nctophonon and magnetoplasmon excitation spectrum of 
these phases. It is therefore conceivable that quantum 
fluctuations, in the form of zero-point energy, may alter 
the cohesive energies and therefore change the structure 
of the phase diagram of the 2DES in higher LLs. This 
point will be investigated in Sec. IIVI where we shall write 
an equation of motion for the electron guiding centers in 
the Wigner and bubble crystals, with the goal of finding 
the normal eigenmodes (and hence the zero-point energy) 
of these structures. An essential ingredient in such a cal- 
culation is the elastic matrix of the Wigner and bubble 
crystals, that we shall determine in Sec. IIIII below. 



where we introduced the effective interaction potential 
U m m' between guiding centers of electrons in states m 
and m', which in real space is given by 

t/w(r)= / p m (q)Vf fF (q)p^(q)e iq - r . (12) 

In the following Subsection, we shall use the above ex- 
pression of the interaction potential to study the normal 
modes of the bubble crystal in higher (n > 2) LLs. 



B. Elastic moduli of the electron crystal 

The derivation of the elastic moduli associated with the 
effective interaction potential in Eq. (|12fl proceeds in a 
standard way as follows. First, we evaluate the elasticity 
matrix <& Q ^(R), which is given byi£ 

*a/s(R) - S R ,o da WRfc) - d a dpU(R), (13) 

k 



III. ELASTIC MATRIX AND ELASTIC 
MODULI OF THE WIGNER AND BUBBLE 
CRYSTALS IN HIGHER LANDAU LEVELS 

Having reviewed the basic phase diagram of the 2D 
electron gas in partially filled LLs, we now would like 
to study the elastic properties of the Wigner and bubble 
crystals. But, in order to be able to do so, we still need to 
derive an effective interaction between the guiding cen- 
ters of the electrons, which in turn will allow us to derive 
the elastic matrix of the WC and BC, and hence find the 
compression and shear moduli of these crystalline struc- 
tures. This will be the object of the following Subsection. 

A. Effective interaction potential between guiding 
centers 

Going back to Eqs. ©-0! if we furthermore write 
for the projected density p(q) the decomposition p(q) — 
E m /5™(q)> with 

. _ f dr\ Vnm (r)\ 2 e-^ 

PmW - e - q Hy iLn{q 2 e/2) > ( W > 

then we can rewrite the cohesive energy Ehf in the form 



with U the total interaction potential between bubbles: 
U(r) = J2 U mm >{v). (14) 

m,m' 

Using the following definition of the direct and inverse 
Fourier transformations (here the q integration is carried 
over the first Brillouin zone of the reciprocal lattice, and 
we remind the reader that A c is the area of the primitive 
unit cell of the bubble crystal): 

$a/3(q) = ^$ a/3 (R)e- jq - R , (15a) 

R 

<& Q/3 (R) = i-^^q)^*, (15b) 

and the identity JHJ, one can easily show that $ Q( g(q) can 
be written in the form 

*a/?(q) = jE [(9a + Qa)(g/3 + Q/3)C/(q+Q) 

c Q 

- Q a QpU(Q)]. (16) 



^ Expanding the second term on the right hand side of the 

Ehf = — J]] [/ mra '(Rj — Rj) + U mm i(0), (11) above equation in q around q = leads to the result 



4 



$a/s(q) = -j-qaqpl U{q) 



E 

Q^O 



U(Q) + (Ql + Q%) 



U'{Q) 



QlQ}( 



QU"{Q) - U'{Q) 



Q#0 



QU"(Q) - U'(Q) 
Q 3 



(17) 



Given that the elastic matrix of a two-dimensional trian- 
gular lattice is of the general form 



$a/3(q) = (Cll - C 66 )q a q/3 + C66q 2 6 a /3, 



(18) 



we see that the compression and shear moduli can be 
extracted from the expression of § a p (q) according to: 



cn(g) 

C66(<?) 



®xx{qx*-) 



§xx(q y y) 



(19a) 
(19b) 



Using the above definitions, and the fact that: 



9 /Ms 



(20) 



we obtain (note that these moduli have units of energy; 
in order to obtain elastic moduli in units of energy per 
unit area, one has to divide by the unit cell area A c ): 



Cll 



C66 



{2itMP) 



2(2nM£ 2 ) 



U(q) + J2 

Q#0 



Q#0 



Qx 



2U'{Q) , q2q2 JQU"(Q)-U'(Q) 



Q 



Q 3 



(21a) 
(21b) 



In the above expressions, the reciprocal lattice vectors Q 
for a triangular lattice are given by 

„ 27r /2n — m „ „\ , 
Q = -(- 7r x + my), (22) 

where the lattice spacing a is given by 
/ 4?rA4" \ 1/2 

"(veO ' (23) 

From Eq. I|21a|) . we see that the long wavelength limit 
(q — > 0) of the compression modulus of the bubble crystal 
is given by 

Cll( ^0)~^)i, ,24) 

and shows the characteristic 1/q plasmon behavior, 
in agreement with the well-known result for the two- 
dimensional classical Wigner crystal in zero magnetic 
fields 

Figure|21shows the variation of the shear moduli of var- 
ious electronic crystals vs. v* in the Landau levels n = 2 
and n = 3. For small values of the partial filling factor 



v* , C66 is an increasing function of v* . cqq reaches a broad 
maximum at some characteristic value u^, and shows a 
decreasing behavior for v* > v* n . At some higher value 
of the partial filling factor v* , the shear modulus turns 
negative, indicating an instability of the electronic crys- 
tal. Similar results have been obtained by Cote et alml, 
using a rather sophisticated generalized random phase 
approximation (GRPA) method. Although the precise 
location of the maxima of the shear modulus curves is 
slightly different in their case, the general qualitative be- 
havior of the shear moduli that we find is in good qual- 
itative agreement with the results of Ref. [2J. It is in 
fact quite remarkable that we were able to reproduce the 
qualitative behavior of the more involved GRPA method 
using a simple wavefunction ansatz. This indicates that 
using single-particle non interacting wavefunctions is a 
good starting point approximation to study the physics 
of electronic lattices in quantum Hall systems in higher 
Landau levels. 

The peculiar shape of the shear moduli curves in Fig. 
121 suggests that there may be a universal scaling law for 
the shear modulus, of the form: 

c 66 (z/) = M a f(v*/M?), (25) 
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FIG. 2: (Color online) Upper panel: shear moduli of the var- 
ious electronic crystals (in units of e 2 /et) vs. v* in the n — 2 
Landau level. The solid line is for the Wigner crystal, the dot- 
ted line is for the 2e BC, and the dashed line is for the 3e BC. 
Lower panel: same thing for the n=3 level. The dash-dotted 
line is for the M — 4 bubble solid. 
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FIG. 4: Magnetophonon and magnetoplasmon dispersion 
curves along the boundary T — J — X — Tof the irreducible ele- 
ment of the first Brillouin zone for the WC at v* = 0.18 (upper 
panel) and of the M = 2 BC at v* = 0.25 (lower panel). Both 
plots are for the n — 2 Landau level. Left axes correspond 
to magnetophonon values and right axes to magnetoplasmon 
values, and are all measured in units of ujq = e 2 jehl. 



with a universal function / and some scaling exponents a 
and (3. The identical way in which cqq increases at small 
v* for different values of M indicates that a ~ j3. Figure 
13 shows that, at least for the first peak, this is indeed the 
case, with exponents a — (3 ~ 0.68 for all three values of 
M, M = 1, 2, and 3. 

In Sec. IVII we will use the results of this Section as 
an input to the variational replica theory of Chitra et 
alH and try to analyze the experimental data of Lewis 
et tt/jli for the microwave conductivity of the two dimen- 



sional electron system in the n = 2 Landau level. For 
the moment, we shall turn our attention to the normal 
modes of the Wigner and bubble crystals in higher Lan- 
dau levels. 



IV. NORMAL MODES AND ZERO POINT 
ENERGY OF THE WIGNER AND BUBBLE 
CRYSTALS 
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v7M™ 

FIG. 3: Universal scaling of the shear modulus for the Wigner 
and bubble crystals in the n = 2 Landau level. Here the 
exponent a ~ 0.68. 



We now turn our attention to the derivation of the 
normal modes of the Wigner and bubble crystals. To 
fix ideas, we shall consider the simplest case of the two- 
electron bubble crystal in the n = 2 Landau level, and 
write for the electron guiding centers r m , (m = 0, 1 be- 
ing the angular momentum quantum number) at a given 
lattice site Ri the decomposition r mi = Rj + u mi , where 
u m i is the displacement of the m-th electron from the 
equilibrium lattice position R^. Expanding the HF en- 
ergy of a distorted Wigner or bubble crystal 

i,m<m / 
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to second order in the small displacements u rn i leads to 
the elastic energy 

E « = \H «ma(Ri)^'(Ri - n/)u m ,p(R 3 ), (27) 



equation of motion for the m-th electron at lattice site 
Ri can be written in the form (m* is the effective mass 
of the electron in the host semiconductor and E a p is the 
totally antisymmetric tensor in two dimensions) 



where the generalized elastic matrix $ has the following 
expression 

$™™ (Rj — Rj) = 6ij5 mm > E d a dpU mm i' (R,f.) 

k,m" 

— d a dpU rnm '(Ri — Rj). (28) 

Note that, by contrast to the previous Section, where 
all the electrons within the same bubble were assumed 
to have the same displacement vector u^, here we con- 
sider the more general case of electrons within the bub- 
ble moving independently from one another. This leads 
to an elastic matrix $™^™ which depends on the inter- 
nal quantum numbers m, ml of electrons within bubbles. 
This last clastic matrix is related to the elastic matrix 
<]>a/3 of the previous section through the equation: 



4> 



a/3 



(29) 



m 



3,771' 

- ^£ a /3^-u m;3 (R l ). (30) 



We shall seek a solution to the above equation of mo- 
tion that represents a wave with angular frequency u> and 
wavevector q, i.e. u ma (Hi) = A ma (q)e^ q Ri_w *\ where 
the A m 's are complex coefficients whose ratios specify the 
relative amplitude and phase of the vibrations of the elec- 
trons within each primitive cell. Substituting the above 
expression into Eq. (|3(J|I results in the following secular 
equation for the Wigner crystal: 



*°°(q) + tww £ 



^ y (q) - 



IU>UJ C 
, ,2 



A 0x 

A y 



0, (31) 



We now are in a position to find the normal modes of the 
bubble crystal. In the presence of a magnetic field, the 



while for the M = 2 bubble solid the secular equation is 
given by 



/ *2°(q) 



*2i(q) 



*Si(q) *£(q) 



$°°(q)+iww c 
V $2},(q) *£(q) 



$°°(q)-iu;a, c §°J(q) 

*Si(q) *ii(q) 

l»oo 



tUJLUc 



$™(q)-^ 2 $°J(q) 

*il(q) - ' 



(32) 



where we defined 

(q) = ^ E (r)e" Jq - Rl , (33) 

i 

and used the fact that $^'(q) = $™™'(q) = ^ m (q). 

The above secular equations have a non-vanishing so- 
lution (for the A ma 's) only if the determinant of the 
secular matrix is zero. This leads to two types of solu- 
tions for the eigenmodes: 23 magnetophonon modes, with 
eigenfrequencies which vanish like q 3 ^ 2 as q — > 0, and 
magnetoplasmon modes which tend to a finite limit as 
q — ► 0. In Fig. 0] we show our solution for the magne- 
tophonon and magnetoplasmon modes of the WC and the 
M = 2 BC in the n = 2 Landau level. As is expected, 
the eigenmode spectrum for the 2e BC has four branches 
(compared to two for the WC). We note that the magne- 
tophonon dispersion curves we obtain are identical to the 



ones obtained by Cote et al& within the time-dependent 
HF approach. To our knowledge, however, the magneto- 
plasmon modes of electronic crystals, shown here in Fig. 

have not been previously studied in higher LLs. 

We now turn our attention to the calculation of the 
zero point energy of the Wigner and bubble crystals. To 
this end, we shall use the approach of Cunningham, 17 
whereby one evaluates the energy of the magnetophonon 
and magnetoplasmon modes at a given number N p of 
predefined points q^ within the irreducible element of the 
first Brillouin zone for the triangular lattice, with appro- 
priate weights on assigned to each one of these points. 
The zero-point energy (per electron) for the M-electron 
BC, in units of e 2 /e£, is then given by 

2M n p 

7 = 1 1=1 
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V* EwC E2eBC 

0.05 0.362825 0.360756 

0.10 0.368391 0.360716 

0.15 0.379527 0.362053 

0.20 0.382453 0.365264 

0.25 0.367643 0.368439 

TABLE I: Zero point energy (per particle) of the WC and 
2e BC for different partial filling factor values v* . Energy is 
measured in units of e 2 /el. 



where luo — e 2 /eM. The resulting zero-point energies for 

.18 



N p = 6 (using the special points and associated weight; 



igtne 

given in Ref. ll7T ) are shown for a range of partial filling 
factor values in Table. [I] We notice that for v* < 0.20 
the zero-point energy of the WC always exceeds the cor- 
responding quantity for the two-electron BC, which sug- 
gests that quantum fluctuations will shift the transitional 
filling factor v\ 2 between the WC and the 2e BC toward 
smaller values. However, this shift in the transitional val- 
ues of v* , of order 0.02, is rather small, and has no sub- 
stantial effect on the overall phase diagram of the 2DES 
in higher LLs. 



V. EFFECT OF FINITE SAMPLE THICKNESS 
AND OF LOWER LANDAU LEVELS 

Let us now investigate the effect of finite sample thick- 
ness and of lower Landau levels on the energetics of the 
Wigner and bubble crystals in higher LLs. To this end, 
we shall replace the bare Coulomb potential in Eqs. I|4a|) - 
(|4b|) with the following effective interaction 



v(q) 



2TTe 2 



-\qi 



(35) 



The parameter A models a finite thickness sample, and is 
generally taken to be of order unity^ The effect of lower 
LLs is encoded in the wavevector-dependent dielectric 
constant e(q), for which we shall use the following ex- 
pression, due to Aleiner and Glazmar 



,20 



e{q)=e{l + ^-[l- J$(qR c )}}, 
I qa B > 



(36) 



with cib = h 2 e/m*e 2 the effective Bohr radius and R c — 
V2n + It 

The form of the effective Coulomb interaction of Eq. 
(|35|l makes it difficult to find an analytic expression for 
the Fock part of the interaction potential, Eq. (|4bjl . To 
evaluate the cohesive energy, the q integration in this last 
equation is performed numerically. In Figs. © and 10, 
we plot the cohesive energies of the Wigner crystal, bub- 
ble phases and stripe states in the second (n = 2) and 
third (n — 3) Landau levels, respectively. In each of these 
two figures, the upper panel corresponds to A = 0, while 
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FIG. 5: (Color online) Cohesive energies of the Wigner and 
bubble crystals in the n — 2 LL in units of e 2 /e£, with the in- 
teraction potential Vhf resulting from the effective Coulomb 
potential in Eq. 1351 . The solid line is for the Wigner crys- 
tal, the dotted line is for the two-electron bubble phase, the 
dash-dotted line is for the M = 3 BC, and the dashed line is 
for the stripe phase. The upper panel is for A = 0, and the 
lower panel is for A = 1. 



the lower panel corresponds to A = 1. The screening 
by lower LLs has the most drastic effect on the cohesive 
energy, shifting it upwards by about 50 percent of its 
bare value, and finite sample thickness also tends to shift 
the cohesive energies up, although in a less pronounced 
fashion. It is noteworthy that these shifts in the cohe- 
sive energy curves do not alter the overall phase diagram. 
As can be seen from Fig. the transitional values v\ 2 
from WC to the 2e bubble phase, and 3 from M = 2 
to M = 3, are slightly shifted downwards by screening 
alone (we find v\ 2 ~ 0.2 and v\ 2 ~ 0.35 with screening, 
while v\ 2 — 0.22 and v\ 3 ~ 0.37 without screening). Fi- 
nite thickness effects on the other hand shift transitional 
filling factors upwards, restoring them to values that are 
very close to their bare values. In all cases, these shifts 
do not affect the overall phase diagram, which remains 
the same as in the ideal case of zero sample thickness and 
no screening by lower Landau levels. 

We have also calculated the shear moduli of the Wigner 
and bubble crystals in presence of screening and taking 
finite thickness effects into account. As we mentioned 
above, in the presence of screening, it becomes difficult 
to obtain an analytic expression for the Fock part of the 
HF interaction potential, and the method of Sec. IIII Bl 
becomes impractical. A more efficient way to extract 
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FIG. 6: (Color online) Same as in Fig.|SJ but for Landau level 
n = 3. The thick dash-dotted line is for the bubble state with 
M = 4. 

the elastic coefficients of that case consists in finding the 
cohesive energy of a distorted crystal for a given uniform 
deformation. The elastic energy will then be given by 
the excess energy of the deformed crystal with respect 
to the reference undistorted state. For two-dimensional 
uniform deformations, such that the displacement vector 
u(r) is given by 

u a = u a ,pxp , (37) 

with constant coefficients u Qj/ g, it can be shown2ithat the 
reciprocal lattice also experiences a homogeneous defor- 
mation, with the deformed reciprocal lattice vectors {Q} 
given in terms of the original ones {Q}, to first order in 
the small displacements {u}, by: 

Qa = Qa ~ Up.aQp. (38) 

To find a given elastic modulus, we calculate the HF en- 
ergy of the corresponding distorted state (with distortion 
amplitude Uo), and extract the elastic constant from the 
excess energy Ehf(uo) — Ehf(0)- For example, to ex- 
tract the shear modulus cqq we calculate the HF energy 
Ehf{uq) of the distorted crystal using the following re- 
ciprocal lattice vectors: 

Qx = Qx - u Q y , (39a) 
Q y = Q y , (39b) 

corresponding to the shear deformation polarized along 
y such that u(r) = UQxy. The shear modulus (in units 
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FIG. 7: (Color online) Upper panel: Shear moduli (in units of 
e 2 /e£) for the WC (solid lines) in the n = 2 LL with the inter- 
action potential Vhf resulting from the effective Coulomb po- 
tential in Eq. 1351 . The succession of curves is for A = 0, 0.5, 1 
and 2. Lower panel: shear moduli of the WC and for the 2 
and 3 electron bubble solids (from bottom to top) in presence 
of screening by lower LLs and with A = 1 . 

of energy /particle) is then given by: 

c 66 = \ [E HF K) - E HF (0)] . (40) 

Note that this procedure allows us to extract the elastic 
moduli without Taylor expanding in the small displace- 
ment vectors {u}, which allows us to better handle anhar- 
monic effects which might be important in case quantum 
fluctuations happen to be larg o 22 ! 2 ? (which might be the 
case with the softer lattices in presence of screening by 
lower LLs). Note also that, since cohesive energy calcula- 
tions involve reciprocal lattice sums that avoid summing 
over the (diverging) contribution of the Q = Hartree 
term, one should add to the results for the compres- 
sion modulus en obtained by the above method a term 
-j-Uiiiq) to obtain the correct compression moduli of the 
electronic solid. 

In order to test the validity of the above procedure, 
we have calculated the shear moduli of the Wigner and 
M = 2 bubble crystals in the n — 2 Landau level. As 
can be seen in Fig. (upper panel), we obtain the same 
shear modulus for the WC as in Sec. IIII Bl where we 
used a Taylor expansion of the interaction energy in the 
small displacements {u(r)}, which is a consistency check 
on our calculations. In Fig. \7\ we also show the results 
we obtain for the shear modulus of the WC in presence 
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FIG. 8: Universal scaling of the shear modulus for the Wigner 
and bubble crystals in the n = 2 Landau level in presence of 
screening and finite thickness effects. Here we used the value 
A = f for the screening strength, and the value a ~ 0.75 for 
the scaling exponent. 



of screening and for various values of the parameter A. 
It is seen that screening by lower LLs alone (i.e., with 
A = 0) considerably softens the shear modulus with re- 
spect to the unscreened case, and higher values of A lead 
to even smaller values of the shear modulus of the WC. 
In the lower panel of Fig. [7\ we plot the shear moduli of 
the WC and of the 2 and 3 electron bubble states in the 
n = 2 LL in presence of screening and using A = 1. As 
it can be seen, the shear moduli of these various types 
of crystals are affected by the screening of the Coulomb 
interaction between electrons in similar but noncquiva- 
lent ways, with the strongest effect impacting the WC. 
This is quite understandable, given that a shear defor- 
mation of a BC does not affect the part of the cohesive 
energy coming from electrons inside the bubbles. It is to 
be noted also that the universal scaling of Eq. 1)25(1 holds 
only approximatively in the screened shown in 

Fig. [5J where we used a — f3 ~ 0.75. We shall discuss 
some of the implications of our results in Sec. IVII 



VI. LINK WITH EXPERIMENT: 
CONSEQUENCE FOR THE MICROWAVE 
CONDUCTIVITY OF TWO DIMENSIONAL 
BUBBLE CRYSTALS 

We now would like to discuss the experimental impli- 
cations of our results. Our investigation of the dynam- 
ics of bubble solids stemmed primarily from a desire to 
understand recent microwave conductivity experiments 
by Lewis et a/^i in the second (n = 2) Landau level, 
where the appearance of a second resonance peak around 
v* ~ 0.16 was attributed to the formation of a bubble 
phase, which coexists with a Wigner crystal over a range 
of values of v* . In what follows, we shall attempt a qual- 
itative analysis of the experimental results of Lewis et 
oZJi in light of the results of the elasticity theory pre- 
sented in this paper combined with the theoretical pre- 
dictions of Refs. Illll 



According to Chitra et al.^ the resonance frequency 
due to pinning of a two dimensional Wigner crystal in a 
strong magnetic field is given by 



(41) 



Pm^c 



where p m = m/(-Ka 2 ) is the mass density, uj c — eB/m 
is the cyclotron frequency, and where £ is given by (we 
neglect an overall numerical constant of order unity): 



(42) 



In the above expression, £ is the greatest of I and the 
correlation length £<j of the disorder potential. On the 
other hand, R a is the length scale at which electron dis- 
placements become of order a, and is given by (here A is 
the variance of the random pinning potential, whose dis- 
tribution was assumed by Chitra et al. to be Gaussian 
and of short range): 



Using Eq. (g3J) into Eq. we obtain: 

A 



and hence 



p m uj c c 66 ^a 2 ' 



(43) 



(44) 



(45) 



Using the fact that p m = m/ira 2 , we see that p m a 2 = m, 
and hence: 



LO p {v*) 



muj c Qc 66 



(46) 



In the above equation, m is the mass of electrons within 
a unit lattice cell, and hence m = Mm* , where m* is the 
effective electron mass in the host semiconductor. Wc 
thus obtain: 



LO p {v*) 



Mm*uj c Qc 66 



(47) 



The above result implies that the v* dependence of the 
resonance frequency lo v arises mainly from the depen- 
dence of the effective shear modulus cqq on the partial 
filling factor v* (the dependence of £o and on v* be- 
ing rather weak in Landau levels of index n > 2). We 
thus see that knowledge of the variation of Cee vs. v* in a 
given Landau level and for a given bubble crystal should 
allow us to easily infer the variation of lu p vs. v* . Con- 
versely, analysis of experimental data of lu p vs. v* using 
the above expression (and our Hartree-Fock results for 
c m( v *)) should allow us to identify Wigner crystal and 
bubble phases, and transitions between these phases. 
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We now turn our attention to the experimental results 
of Ref. 0, and more specifically to Fig. 3 in this last 
reference, where the authors plot experimental resonance 
frequencies of the two-dimensional electron system in the 
ri = 2 Landau level as a function of the filling factor v. 
At small v* , there is only one resonance peak, whose 
frequency decreases with increasing v* , with a minimum 
around v* ~ 0.19. The frequency of the above peak 
then slightly increases, until the peak disappears around 
v* ~ 0.27. This behavior is in qualitative agreement 
with the behavior of a Wigner crystal if we use our shear 
modulus of Fig. [21 (Wigner crystal in the n = 2 LL) in 
conjunction with Eq. (|47|l . 

Figure 3 of Lewis et al. also shows the peak frequency 
of a second resonance peak which appears at v* > 0.16. 
The v* dependence of the frequency of this second res- 
onance peak, however, does not seem to be fully consis- 
tent with Eq. I|47|l above and with the shear moduli of 
the M — 2 bubble phase plotted in Fig. [3 Indeed, in the 
range 0.23 < v* < 0.35, from the decreasing behavior 
of the shear modulus of the bubble phase in Fig. [21 we 
expect an increasing behavior of the peak frequency u) p 
vs. v* . The experimental result of Ref. ^] show a de- 
creasing peak frequency in this range of filling factors, in 
disagreement with the theoretical considerations above. 
Equally important is the fact that Eq. (|47|l predicts a 
factor of 1/M difference between the frequencies of dif- 
ferent bubble states, while experimentally it seems that 
a factor 1/y/M is observed instead. 

It is to be noted that, in their recent replica study of 
bubble pha ses pinned by random disorder, the authors 
of Ref. [24] claim to have been able to successfully de- 
scribe the v* dependence of the resonance peak of both 
Wigner and bubble crystals, except on the region be- 
tween v* — 0.16 and v* — 0.28 where the two phases 
are assumed to coexist. We therefore may speculate that 
the disagreement we find in the present study may be 
due to a possible breakdown of Eq. (|4ip. which in fact 
is derived by Chitra et al. 12 by linearizing nonlinear self- 
consistent replica calculation, and that a careful and ex- 
haustive treatment such as the one in Ref. 0, where the 
full nonlinear replica equations are numerically solved in 
a self-consistent way, can lead to better agreement be- 
tween theoretical predictions and experimental observa- 
tions. 

We now briefly comment on the effect of screening by 
lower LLs and of finite sample thickness (as modeled by 
the parameter A of the Zhang-Das Sarma potential of the 
previous Section) on the microwave response. Since we 
have shown that both effects tend to reduce the shear 
modulus with respect to the bare case, we conclude that 
the resonance peak will tend to shift toward higher fre- 
quencies as the sample thickness or Landau level index 
is increased, an effect that may be testable experimen- 
tally. More importantly, in presence of screening and 
for finite thickness samples, we find that there is a finite 
gap between the shear moduli of the Wigner and bubble 
crystals, and that gap might accentuate the finite gap 



observed in the coexistence region between the pinning 
frequencies of the Wigner and the 2e bubble crystals, 
making it more pronouncedii In fact, it would be inter- 
esting to do a full self-consistent calculation, like the one 
done in ref. 12 ll in presence of screening and finite thick- 
ness effects, and see whether the experimentally observed 
gap between the pinning frequencies in the coexistence 
region can be reproduced theoretically. 



VII. CONCLUSIONS 

In conclusion, in this paper we have calculated the 
cohesive energies, shear moduli, and normal (magne- 
tophonon and magnetoplasmon) modes of the Wigner 
and bubble crystals in higher LLs. Going beyond previ- 
ous treatments, we have studied the effects of screening 
by lower Landau levels as well as the case of finite sample 
thickness. We found that both effects reduce the cohe- 
sive energies and elastic moduli, with the screening by 
lower Landau levels having the most pronounced effect. 
The transitional values of the filling factor as well as the 
overall nature of the phase diagram remain, however, un- 
changed by both effects. We have also examined the elec- 
tromagnetic response of the Wigner and bubble crystals, 
and have noticed that, while the v* dependence of the 
shear modulus of the WC, when used in standard theo- 
ries of the electromagnetic response of electronic crystals, 
is in good qualitative agreement with experimental data, 
the shear modulus of the 2e BC has a v* dependence 
which, in contrast, is in disagreement with the experi- 
mentally observed behavior. This led us to suggest that 
a full numerical solution of the non-linear self-consistent 
replica equations^ may be needed in order to understand 
the electromagnetic response of bubble crystals, by con- 
trast to Wigner crystals for which a linearized estimate 12 
of the resonance frequency may be adequate. 

Before closing, we briefly comment on an issue that has 
not been previously addressed in the literature, and that 
is the issue of gauge invariance. In principle, all physical 
observables are gauge-independent. However, HF being 
an approximate (i.e. not an exact) theory, different HF 
derivations done in different gauges will not necessarily 
give the same results. As a check to our calculations of 
the cohesive energies of the Wigner and bubble crystals, 
and to make sure that the same qualitative phase dia- 
gram that we obtained above in the Landau gauge is not 
affected by a different choice of gauge, we have rederived 
the HF interaction potentials and the phase diagram of 
the 2D electron system in the symmetric gauge, using 
a HF approach similar to the one used a long time ago 
by Maki and Zotosi^ to study the properties of the two- 
dimensional Wigner crystal in the lowest Landau level. 
In Appendix [21 we derive expressions for the interac- 
tion potential U mm > (r) between electrons of angular mo- 
menta m and ml ', with m,m' = 0,1,2 by finding the 
average Coulomb energy between symmetric-gauge elec- 
tronic wavefunctions in the n — 2 Landau level. These 
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interaction potentials have the same qualitative behav- 
ior as the potentials derived within the Landau gauge 
HF approach of the text, and lead to very similar cohe- 
sive energies as the ones predicted above, although with 
slightly different numerics. 
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densities and Hartree-Fock potentials that we used in 
deriving our cohesive energies. From the usual expres- 
sions of the noninteracting wavefunctions vs„ m (r) (see 
Eq. i|C2(l . and the definitions we obtain that the 

projected densities po(q) and pi(q) are given, both in 
n = 2 and n = 3, by (throughout this Appendix, q stands 
for the dimensionless quantity ql ) 



p (q) - e-f\ (Ala) 
pM) = (l-lq 2 )e-« 2 /\ (Alb) 



APPENDIX A: PROJECTED DENSITIES AND 
HF POTENTIALS 

In order for our paper to be self-contained, in this Ap- 
pendix we give the explicit expressions of the projected 



On the other hand, performing the integrals in Eqs. 14a|) - 
(|4b|) , we obtain the following Hartree and Fock potentials 
in LL n = 2: 



(2-Ke 2 l\ e~ q2 / 2 , 9 4n2 



(82 - 52q 2 + Uq 4 - I0q 6 + q 8 )l {^-) ~ c? 4 (30 



e / 128 

while for LL n — 3 the Hartree and Fock potentials are given by: 



V, 



n(q) = ( 
V>(q) = - 



27re 2 A e-« 2 / 2 



e / 2304 9 

4608 V 2 



(48 - 72q 2 + 18g 4 - q 6 



g 4 (l539 - 828<? 2 + 224 9 4 - 24q 6 + g 8 )li(^-) 



(2646 - 2430g 2 + 2889 9 4 - 1236g 6 + 270g 8 - 26g 10 + g 12 )l (^-) 



(A2a) 
(A2b) 

(A3a) 
(A3b) 



r 



APPENDIX B: COHESIVE ENERGY OF THE 
STRIPE PHASE 

For completeness, in this Appendix, we briefly re- 
view the cohesive energy of the stripe phase of the two- 
dimensional electron system. Following Fogler et a/i and 
Goerbig et alA, we write the cohesive energy per particle 
E co h m the form 



£c 0ft =U$>^(g)|A(g)f 



(Bl) 



where ns — 1/(2tt£ 2 ) is the flux density, and where 
A(q) = p(q)/(ngA) is the Fourier transform of the lo- 
cal guiding center filling factor v* (r) (A is the area of the 
sample). For a stripe state, v*{r) = JY Q(\x — X{\ —a/2), 



with a the width of a given stripe, and Xi the location 
of the center of the i-th stripe. If we denote by L s the 
stripe periodicity (such that X{ — iL s ), and note thali 1 ^ 
v* = a/L s , then we obtain after a few manipulations that 
the cohesive energy per particle of Eq. I|B1|1 is given by 



n B 
2ttV 



HF 




sin 2 \ixv* n) 



(B2) 



where the prime on the sum sign indicates that the (di- 
verging) Hartree contribution for n = is excluded from 
the summation (the n — Fock term being included). 
The above expression is then minimized with respect to 
L s , with the optimal stripe periodicity found^ to be ap- 
proximately given by L 
L 



2.74£V2n + 1 in n 



t 2.1U^J2n + 1 in n = 2 and 
3. The cohesive energy in the 
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FIG. 9: Upper panel: Interaction potential Uoi(r) (in units of 
e 2 /£) vs. r (in units of £) as obtained from the HF approach 
in the symmetric gauge, Eqs. l(J12el - ^C12fl> . in the n — 2 
Landau level. The dotted line corresponds to the numerator 
iVoi(r) alone. Lower panel: Effective interaction potential 
Uoi(r) in the Landau gauge. 



stripe phase is then given by the optimal value E co h{L*). 




FIG. 10: Interaction potential Uoo (in units of e 2 /£) vs. r (in 
units of £) in the n = 2 Landau level. The solid line is the 
Landau gauge result, while the dashed line is the result for 
Noo{r) in the symmetric gauge, Eq. (IC12al> . The dotted line is 
the whole function Uoo(r) — Noo(r)/Doo(r). (Note that this 
last function has no physical meaning at small separations, i. e. 
for r < £, due to the Pauli exclusion principle for Fermions.) 



where <£> nm R(r) is the symmetric gauge wavefunction 
with angular momentum m which is located around the 
lattice site R fc , (p nm n.(r) = ip nm (r - R), with 



CO (r) - C (W n - ml P ^-m)(> p -r 2 /±e 2 



x L 



| n — m | 

(n+m— \n— m\)/2 \ 2£2 



(0^2)' 



(C2) 



where we denote by 9 the polar angle of vector r, i.e. 
r = (r cos 9, r sin 9), and where the normalization con- 
stant C nm is given by: 



APPENDIX C: GAUGE DEPENDENCE OF THE 
COHESIVE ENERGIES: HARTREE-FOCK 
APPROACH IN THE SYMMETRIC GAUGE 

In this Appendix, we want to confirm the qualitative 
features of the phase diagram of electrons in higher LLs 
derived in the text using a different gauge for the ap- 
plied magnetic field, namely the symmetric gauge A = 
y(— yx + xy). To this end, we shall use a HF ap- 
proach which is similar in spirit to the approach of Maki 
and Zotosi^ to derive the effective interaction potentials 
U mm i{r) in the second (n = 2) Landau level. Let us con- 
sider two electrons with angular momenta mi and 7712 
and with guiding centers located at lattice sites Ri and 
R2, respectively. We shall use for the wave function of 
the two-bubble system the following antisymmetric com- 
bination 

*(ri,r 2 ) = -y=((Pnm 1 R.i{ri)fnm 2 R. 2 { r 2) 

~ Vnm 1 R 1 (r2)V ; 'nm 2 R 2 (ri)), (CI) 



On 



2irm\ 



27rn ! 



if m > n, 



if n > m. 



(C3) 



The average potential energy for this state is given by 
(we here for simplicity set e = 1): 



,(Ri, R 2 ) — 



ri— r 2 | 



I*) 



(*|*) 



(C4) 



Note that the above expression is valid for arbitrary Ri 
and R2 only if mi ^ mi. In the special case mi = m-i-, 
because of the Pauli exclusion principle for Fermions the 
above expression for U mini2 only makes sense at separa- 
tions |Ri — R2I larger than the typical extension of the 
wavefunctions tp nm (r), which is of order £. 

Now, the numerator of the above expression is given by 
(to simplify the notation, in the rest of this Appendix we 
shall drop the LL index n from the ip nm wavefunctions): 



13 



7V mi , m2 (R 1; R 2 )=e 2 / | _ | [I Vmi,Ri(ri)| 2 |<^m2,R2( r 2)| 2 ~ <Pmi Ri( r lVm 2 R 2 (l"2)^m 2 ,R 2 {Tl)<p mi ,Ki ( r 2)] ■ (C5) 
J ri — r 2 1 



The denominator on the other hand is given by 



£> mi , m2 (Rl,R 2 ) = 1 - / *l*2^m 1 3l( r l)V'm 1 ,R 1 (l"2)<y3 m2: R2( r 2)^m 2 ,R 2 (ri) =1 ~\ J ^ VW.Ri (*1 )</4 2 ,R,- ( r l) | 



Now, using the fact that 



(C6) 



^a(r) = fe(r - R) = [ ^ m (q)e- iq - (r - R) , 



(C7) 



one can show, after a few manipulations, that the Hartree part of N mim2 is given by 



< im3 (R i ,R J .) = ^^( | dr 1 | ¥ , mi (r 1 )| 2 e 4q - ri )( / dr 2 |^ m2 (r 2 )| 2 e -^) 



2 p -iqT2^giq-(Ri-Rj)_ 



(C8) 



Similarly, the Fock part of N mirrl2 can be written in the form 



N miTO2 (Ri,Rj) - - 



27re 2 



(r)^ m2 (r + Rj - R^e' 1 



(C9) 



r 




Eqs. HC8 MC10|1 show that both Hartree and Fock 
parts of N mim2 as well as the denominator D mirrl2 de- 
pend only on the difference R; — Rj and are thus trans- 
lationally invariant as they should be. Using for the 
Vm.R's the wavefunctions for noninteracting electrons in 
the n-th LL, Eq. I|C2|) . and performing the integrations in 
Eqs. HC8fl - l|C10|) (which is most efficiently done in carte- 
sian coordinates), we obtain that the above potential en- 
ergy, Eq. I|C4(1 . can be written in the form: 



FIG. 11: Interaction potential Un (in units of e 2 /£) vs. r (in 
units of €) in the n — 2 Landau level. The solid line is the 
Landau gauge result, while the dotted line is the symmetric 
gauge result for Nu(r), Eq. lC12ct . 



U mi m 2 (Rl — R2) 



e 2 N mim2 (Ri — R 2 ) 
( D mim2 (Ri — R 2 ) 



(Cll) 



Finally, the denominator D mim2 can also be written in 
the form 

D mima (Ri,Rj) = 1 - I J drip* mi (r)(f m2 (r + R^ - Rj)| 2 . 

(CIO) 



where the functions N mim2 and D mim2 are given, for 
n = 2 and mi,m 2 = 0,1, by (Io and Ii are modified 
Bessel functions) 
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N Q0 (r) 



/7re 4 
32768 



(9360 - 2944r 2 + 560r 4 - 32r 6 + r 8 ) 



65536 



(18720 + 992r 2 + 240r 4 + 12r 



)I (-) - r 2 (2432 + 312r 2 + 16r 4 + r 6 )!^-) 



A»(r) = l-^(l28-32r 



2 , ^ 



(C12a) 
(C12b) 



Nn(r) 



y 7re 4 
2097152 



(520384 - 400512r 2 + 150672r 4 - 24448r 6 + 1964r 8 - 72r 10 + r 12 ) 



'ire 



1040768 + 27840r 2 + 47120r 4 - 6656r 6 + 904r 8 - 44r 10 + r 12 )l (— ) 

8 

2 



4194304 

r 2 (169088 + 35760r 2 - 3904r 4 + 752r 6 - 40r 8 + r 10 )li(y ) 

r 2 



On (r) = 1 - 
N i(r) = 



1048576 

ire 5" 



262144 



( - 1024 - 384r 2 - 40r 4 + r 6 ) , 

(10784 + 16144r 2 - 5152r 4 + 816r 6 - 46r 8 + r 10 ) 



(C12c) 
(C12d) 



1 



(128192 + 16976r 2 + 464r 4 + 268r 6 - 12r 8 + r 10 )I (— ) 



- r 2 (24976 + 1440r 2 + 244r 4 - 8r 6 + r 8 )Ii (— ) 



D m (r) = 1 



262144 



(128 - 32r 2 



(C12e) 
(C12f) 



r 



A plot of the function Uoi(r) (see Fig. EJl shows the 
same qualitative behavior as the one obtained within the 
Landau gauge, Eq. 1)12(1 . with, namely, a negative curva- 
ture at the origin and a local minimum at a finite value 
r = 7q ~ \A\l. Because the wavefunctions (fio(r) an d 
¥>i(r) are orthonormal, the denominator Dqi{t) is al- 
ways very close to unity, and the potential Uq\{t) is very 
well approximated by the numerator iVoi(r). The situa- 
tion is, however, markedly different for the pairs of wave- 
functions (fio(r — R) and <fio(r — R'), and </?i(r — R) and 
i^i (r — R'), which have very strong overlap at small sep- 
arations |R — R'|, with the consequence that the denom- 
inators -Doo( r ) and Dn(r) become very small for small r, 
leading to nonsensical results for the interaction poten- 
tials U 00 (r) and t/n(r) (see Fig. EU- As it can be seen 
in Figs. ifTUl ITT|) . we find that a good approximation to 
the true Hartree-Fock interaction potentials Uoo(r) and 
U\\{r) is obtained by keeping the numerator contribu- 
tions only. Similar results are obtained for the interac- 
tion potentials Uo2(r), Un{r), U^r), but we choose not 
to write down the explicit expressions of these potentials 
here for brevity. 

Using the above effective interaction potentials be- 
tween guiding centers leads to the cohesive energies for 
the Wigner and M = 1 and 2 bubble crystals in LL n = 2 
shown in Fig. 1121 As it can be seen, the qualitative pic- 



ture of the relative behavior of these cohesive energies is 
the same as in the approach of the text, which is based 
on the Landau gauge. We thus see that the phase dia- 
gram of electrons in higher LLs is quite robust, and that 
a particular choice of gauge is likely not to have an ef- 
fect on the qualitative results of HF calculations, even 
though slight quantitative discrepancies between calcu- 
lations carried out in different gauges might exist. 



APPENDIX D: NUMERICAL VALUES OF 
PHYSICAL QUANTITIES 

In this Appendix, we give a summary of numerical val- 
ues of physical quantities we use to calculate the normal 
modes in Sec. IIVI Throughout this Appendix, we shall 
use CGS units, where the speed of light in vacuum c, 
Planck's constant h, the electron charge e and effective 
mass in the host semiconductor m* are given by 



= 3 x 10 10 cm/s, 

= 6.6262 x 10~ 27 erg - s 

= 4.80 x 10~ 10 esu, 

= 0.067 x 9.11 x 10~ 28 

= 6.103 x 10~ 29 g. 



(Dla) 
(Dlb) 
(Die) 

(Did) 



15 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



FIG. 12: (Color online) Cohesive energy (in units of e' 2 /£) vs. 
v* for the Wigner crystal (solid line), for the 2e bubble solid 
(dotted line), and for the 3 electron bubble crystal (dashed 
line) in the n — 2 Landau level, as obtained from the wave- 
function approach of this Appendix, which is based on the 
symmetric gauge. The relative behavior of these cohesive en- 
ergies is very similar to the behavior derived in the approach 
of the text based on the Landau gauge. 



Following Ref. |3|, for the dielectric constant of the host 
semiconductor, we use the value e = 12.9, relevant to 
GaAs, and take for the electron density the typical value 
no = 3.2 x 10 11 e _ /cm 2 . The resulting expression of the 
magnetic field B = nohc/ev in terms of the filling factor 
v is as follows 



B 



13.24 



x 10 Gauss. 



(D2) 



The magnetic length £ = y/hc/eB is then 

1 = 0.7052 x lO- 6 y/U cm, (D3) 

while the effective Bohr radius as — h 2 e/m*e 2 is given 
by 



a B = 1.0188 x 10~ b cm, 
so that the ratio as/i is given by 
a B _ 1.4426 



(D4) 



(D5) 



Now, the electrostatic energy scale e 2 /e£ is such tha^2a 
e 2 24046.7 



Het 



GHz. 



(D6) 



Using Eq. (|D2(1 we obtain that the cyclotron frequency 
lu c = eB/m*c is given by 



3474L2 
lo c = GHz, 



(D7) 



so that the ratio of the electrostatic to magnetic energy 
scales is given by 



2 /e£) 



= 0.692 y/u. 



(D8) 
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